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OF 
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BY 
SLIDING MODE 


MLR. JAFART, M.R. ZARRABI'”, S. EFFATI? 


Abstract. Jn this paper a pray-predator system that called self-remission tumor is 
considered, and a new approach in order to stabilizing the unstable equilibrium 
points of self-remission tumor system with sliding mode control is introduced. 


The stability analysis of the biologically feasible equilibrium points is presented by 
using the Lyapunov function. 


A Lyapunov function is supposed for designing a sliding surface (SS). 


Lyapunov function is constructed to establish the global asymptotic stability of the 
uninfected and infected steady states by describing sliding surface (SS), after that 
by considering the derivation of SS as zero, someone can achieve the equivalent 
control that inbreed system stays on SS and tends to equilibrium point in infinite 
horizon. 


In addition, numerical examples are presented to illustrate the effectiveness of the 
proposed method. 
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1. Introduction 


Cancer is one of the greatest killers in the world and the control of tumor growth 
requires special attention [9]. 


The mathematical modeling of cancer self-remission and tumor has been 
approached by a few numbers of researchers under using a variety of models over 
the past decades [8, 9, 13, 18]. 


Many authors have discussed the problem of the chaotic behavior and stability 
analysis of some biological models such as cancer and tumor model, genital 
herpes epidemic, stochastic lattice gas prey-predator modes [7, 8, 13] and many 
other models. 
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Optimal treatment scheduling of prey-predator system using a control theoretic ap- 
proach is the subject of substantial research activity. In [1,2,6,14,16], open-loop type of 
optimal controllers is designed by using the Pontryagins Maximum Principle (PMC). A 
major drawback of open-loop optimal controllers is their lack of robustness against dis- 
turbances or model uncertainties. In fact, tumor dynamics are poorly known, this leads 
to model inaccuracies and parameter uncertainties. Also, another source of disturbances 
may arise from immune system fluctuating or immune effect of a co-infection, in addition 
to the measurements errors and estimation errors when using an observer to estimate the 
unmeasured states. Therefore, the design of optimal treatment schedules based on open- 
loop optimal controller, may lead to undesired results. To overcome this problem, we 
have to design a feedback controller, which inherits certain robustness to disturbances. 

Feedback control for prey-predator has been studied by [3-5]. A common method for 
using feedback control is sliding mode control that introduced by Emelyanov [10] in the 
late 1960s based on the conceptions of variable structure control (VSC). However, Figge- 
Lotz [11] was the first to present the concept of sliding motion, and Filippov [11] was 
the first to consider the solution of differential equations with a discontinuous right-hand 
side. Filippovs pioneering work still serves as the basis for work in sliding mode control 
which was essentially developed by Utkin [11,12,19,20] and Emelyanov [10], Draenovie 
and their co-workers. 

The present paper is concerned with the problem of design equivalent control of the un- 
stable equilibrium points of cancer self-remission and tumor system. Mathematical model 
of spontaneous tumor regression and progression as a deterministic prey-predator model 
have been constructed. Chaos and stability analysis of equilibrium points of the system 
and their biologically existence conditions are perused and inputs control that stabilizes 
these unstable equilibrium points by using sliding mode control are derived. Eventually, 
the results of this paper complement and extend recently published results by El-Gohary 
[9] and Sarkar [18]. 

This paper is organized as follows: Section 2 briefly introduces self remission tumor 
and its mathematical model. Section 3 addresses the equilibrium points and their stability 
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analysis. In Section 4, equivalent control is designed, and examples are presented in sec- 
tion 5 to illustrate the procedure and its validity of the proposed control design. Finally, 
conclusions are given in Section 6. 


2 The mathematical model 


In this section we start to construct the spontaneous tumor regression and progression 
system using a prey-predator system [9]. We consider the predator is T-lymphocytes and 
cytotoxic macrophages/ natural killer cells of immune system, which attacks, destroys or 
ingests the tumor cell. The prey is the tumor cells which are attacked and destroyed by the 
immune cells. The predator has two states, hunting and resting, and destroys the prey. The 
tumor cells are caught by macrophages which can be found in all tissues in the body and 
circulate round in the blood system [9,18]. Macrophages absorb tumor cells, eat them and 
Telease series of cytokines which activates the resting T-lymphocytes that coordinates the 
counter attack. The resting predator cells can also be directly simulated to interact with 
antigens. These resting cells cannot kill tumor cells, but they are converted to a special 
type of T-lymphocyte cells called natural killer or hunting cells and begin to multiply 
and release other cytokines that moreover simulate more resting cells. This conversion 
between hunting and resting cells result in a degradation of the resting cells undergoing 
natural growth and an activation of hunting cells. To introduce the mathematical model 
we assume that the tumor cells are being destroyed at a rate proportional to the tumor 
cells densities according to the law of mass action. Next we also assume that the resting 
predator cells are converted to the hunting cells either by direct contact with them or by 
contact with a fast diffusing substance produced by hunting cells. We consider that once a 
cell has converted, it will never return to the resting stage and active cells die at a constant 
probability per unit time. Finally we assume that all of resting predator and tumor cells 
are nutrient rich under going mitosis and the tumor cells have a proliferative advantage 
over the normal cells [9]. If M(t), N(t), and Z(t) denote thedensities of tumor cells, 
hunting predator cells and resting predator cells at time t, respectively, the resulting system 
dynamical system can be described the following set of nonlinear ordinary differential 
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equations: 
at =q4+rM(1—k;'M) —aMN, (l) 
dN 
< =BNZ— dN, 
dZ a 


where g is the conversion of normal cells to malignant ones, r is the growth rate of the 
tumor cells, k; is the maximum carrying capacity of tumor cells, ky is the maximum 
carrying capacity of resting cells (also, ky > ka), 8 is the conversion rate of resting cells to 
hunting cells, d; is the natural death of hunting cells, dy is the natural death rate of resting 


cells, s is the growth rate of resting predator cells, a is the rate of predation/destruction 
of tumor cells by the hunting cells. Cancer self-remission and tumor system have to be 
analyzed with the initial positivity conditions M(0) > 0, N(0) > 0, Z(0) >0 [1,2]. 


2.1 Parameters reduction 
We reduce the system parameters (1) by using new variables and new parameters. Now 
consider the following new variables: 
T= Qky't, 1 =ky'M, t2=akq"'N, 13 =ky'Z. (2) 
So we have: 
T=@k;"t => thurs S = ko" 
Now consider the following constants a,i = 1,--- ,6 


a;3=rkg a= Bkykgq-' a3 = kyd,q"* (3) 
ag=skq* as=far* — ag = dakyq* 


Substituting (2) and (3) into (1) we get the following system [9]: 
Z = 14 4,7;(1 — 2%) — 272, (4) 
Ly = QgToT3 — A3To, 
Ty = O4T3(1 — T3) — GgTa7T3 — QgTs, 
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Fig. |. The densities of tumor cells, hunting predator cells, resting predator cells for the parameters 
a, = 0.4, a, = 5.9, ag = 0.1, ag = 0.5, ag = 0.06 and a, = 0.05 and the imitial densities are 
7,(0) = 0.3, 9(0) = 1.5 and z4(0) = 0.5. 


3 Equilibrium points, stability analysis 


In this section will be discussed about all of the positive equilibrium points of nonlinear 


system (4) and analyse their local stability by Lyapunov linearization approach. 
The positive equilibrium points of the nonlinear system 4 are as follow: 


1 4 
Ey — 5a af 1+ am (5) 
At = as 
Fa=[5(+ 1+ ogo 
Ba = {5 -1(01 — 22) + Van 04)” Fay), au(as — a9) — o30e]/030, =}, 
The positive equilibrium points are biologically feasible. Therefore Fo is feasible if ag < 
ay. This means that £5 exists and biologically admissible if the natural death rate of 
Testing cells is less than the growth rate of the resting predator cells. Fy is feasible if 


a3 4% 
[aa(@a — a3) — a0¢]/an05 >O—> 2 + 2 <2 (6) 
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Therefore the equilibrium points E,, E,, FE, under the following conditions are positive 
and biologically admissible [9]: 


@) ay 


3.1 The local stability analysis the equilibrium points 

To investigate the local stability positive equilibrium points under the conditions (7) by 
Lyapunov linearization approach, we obtain the Jacobian matrix J and its eigenvalues at 
these points. 
The Jacobian matrix of the system (4) at the equilibrium state F; is given by 


—/faf+4a, —(14+,/1+ 4)/2 0 
J\z, = 0 a3 0 . (8) 
0 0 Qy — Og 


The eigenvalues of the Jacobian matrix J; are given by 


a -a + 4a;, (9) 


Aya = —43, 
Ais = G4 — Gg, 


From the condition of the (7) we observe that the eigenvalue 4,5 is positive and hence the 

F, equilibrium point is absolutely unstable [15,17]. 

The Jacobian matrix of the system (4) at the equilibrium state F» is given by 
—/faf+4a, —(1+,/1+4)/2 0 


Jie, = 0 (aya —agaq—ay0¢)/ag 0 |- (10) 
0 —@x (4 — Og) / 04 ay — A 


The eigenvalues of the Jacobian matrix Jo are given by 


An = —1/a? + 4a, (11) 
Ang = (@204 — @304 — O06) / a4, 
Agg = Gg — Og, 
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From the condition of the (7) we observe that the eigenvalue Aq is positive and hence the 
F, equilibrium point is absolutely unstable [15,17]. 
The Jacobian matrix of the system (4) at the equilibrium state F3 is given by 
%—k (a; —%2+k)/2a, 0 
Jz, = 0 0 ri: 


a7T2 (12) 
0 —(@gax)/a2 —(@3a4)/a9 


where k = \/(a1 — Za)" + 4a: and Z = |as(a2 — a3) — a206]/a205. The eigenvalues of 
the Jacobian matrix J3 are given by 


Asi = 7 —k, (13) 


Aso = [—@ga@4 + aja} — 40303079] /209, 


Ass = [—agaq — yf 030 — 403030,79]/2a9, 


From the condition of the (7) we observe that the eigenvalues A3;, Ago and Ags have 
negative real parts therefore the equilibrium point F, is asymptotic stable [16,17]. 
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Fig. 2. The cancer self-remission and tumor system converges to stable inside limit cycles for the 
parameters a, = 2.5,a, = 4.5,a, = 0.6,a, = 3.5, a, = 2 and ag = 0.1 and the imitial densities are 
2, (0) = 0.5, 29(0) = 1,273(0) = 1.5 
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Fig. 3. The cancer self remission and tumor three dimensional phase plot which represent three different 
attractors for the values of the system parameters and initial densities 
@, = 5,4, = 15,4, = 0.5,a, = 6.5,a5 = 0.1 and ag = 0.21 and the imitial densities are 
7,(0) = 0.5, 2(0) = 2, 23(0) = 1. 


4 Designing equivalent control 


In this section we will study the problem of design equivalent control of the cancer 
self-remission and tumor model about its equilibrium points. For this purpose of in this 
system as both Z and N cells are from T lymphocyte, we add the control u to system (4) 
in order to reinforce the immune system of body against the tumor cells. Then we will 
use sliding mode approach. Therefore let us consider the coupled system of the cancer 
self-remission and tumor has the form 


Zt = 1+ a,;2,(1 — 7) — 2429, (14) 
To = QgT9T3 — AgTg + Ayu, 
£3 = Q4T3(1 — 13) — GgToT3 — gTs + Qu, 


To study the problem of design equivalent control of the system (14) about the equi- 
librium points £, and Fy we will obtained the equations of perturbed states about these 
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equilibrium points by introducing the following new variables 
m=%1—-T1, M=—T2—T, hs = 73 —Ts. (15) 

where T,i = 1,2,3 denote the coordinates of the equilibrium points FE; and Fy. Sub- 
stituting (15) into (14) taking into account the identities that satisfy by the equilibrium 
points, we get the following system: 

th = (a1 — 2a; 7 — Ta)m — aif — Tim — THM, (16) 

The = (2%3 — 3) + gos + O72 + OU, 

tis = —OsT aM + (04 — Og — O57 y — 2agTs)Ns — O47 — OxtaTs +024, 
In this phase we delete the control agu than third equation. 


u=, Y=, Ms = Ms — (2). (17) 
Therefore we get the following system: 
1 = (a; — 20,7; — Toy — ary? — Try — yx%, (18) 


t= (F)ast + Gg73 — a3)0 + (aaa + Q9To¥3 + G2Y2¥3 + OU, 
ds = Cyyg + Coy2 + Coys + Cayo — a4y3, 


G = (> )las +hg—-o (Pa, + g)F_ — (a2 + 204 + (= )as)a), (19) 
__(@ a 
y= —(=*)" lan +04 + (Jar, 
Cy = [04 ~ 06 — (()ag + a5)%2 ~ 2047), 
a at 
C= ~(G Mee + 2a4 + (ol: 


Consider V (ys, yo, vs) = $(u7 4-5 +45) as the candidate Lyapunov function of the system 
(18). The time derivation of V along the trajectories of the system (18) is 
V = win + vote + ysts (20) 
= (ay — 2ayT; — To)yz — ayy} — Tyy1 Yo — yive 
+(()aat» + a3% — a3)y3 + (S)aaug + QgTa¥o¥s + 0243Us + YooyU 
+Cyyous + Cougs + Caus + Caves — a4y3- 
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First by factorization from linear coefficients of the y and then by factorization from 
nonlinear coefficients of the y, V convert to the following equation: 


V = wo(—Tins — uf + Diys + Doug) (21) 
+to( Dayo + Davos + Dyy3 + au) 
+f (v1, 4s). 
where 
f (us, vs) = (a1 — ays — Ta) uz — (= )*(ag +04 + (S)as)u8 — O4u3, (22) 


Di= (Sas $04 — ag — (Jan n (F yan + as) — (a2 + 204+ ([)as)z), 
Dy= (Naa +2044 (sas); 

Dy = ag — a5 + (asta, 

Dy=—04 - (S)as, 

Ds = ()as, 


The sliding surface can be described as: 
yo = —K(-Ziv — vi + Diys + Days). (23) 
where K > 0. To get the equivalent control we can equal the nonlinear coefficient of the 
y to zero. 
Ueq = —0;'(Dsy2 + Dayays + Dsy3). (24) 
So V can be written as: 
V =—Kyg + f(ui,4s). (25) 
Now, we should assign K in a way that V < 0, therefore the following inequality is 
obtained. For its numerical calculate, we use matlab software. 
K > uy? f (us, us). (26) 


By substituting sliding surface (23) in system (18), the dimension of this system decrease 
into two, and system trajectories approach to the origin by sliding surface (23). So system 
(16) is asymptotically stable in Lyapunov sense. 
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5 Analysis and numerical simulation 


In this section the problem of the numerical solution of controlled nonlinear cancer 
self-remission and tumor system is considered, to exhibit the control of this system by 
using SMC. Numerical examples for controlled cancer self-remission and tumor system 
were carried out for various parameters values and different initial densities. The fol- 
lowing figures display the stabilized behavior of the cancer self remission system about 
the unstable equilibrium points and its equivalent control input (See Figs. 4 and 5). We 
conclude that all of the unstable equilibrium points of the cancer self-remission and tu- 
mor system can be optimally asymptotically stabilized with nonlinear equivalent control 
inputs 04 Ue, and OoU,, under the conditions a, > ag and 3 + % < 2. 

In Figure. 4, 7;, r2 and 1; the densities of tumor cells, predator hunting cells and resting 
predator cells respectively and equivalent control u,, in equilibrium point FE), are dis- 
played against time for the system parameters and initial densities ay = 0.4, ag = 5.9, 
a3 = 0.1, ag = 0.5, ag = 0.6, ag = 0.05, ay = 0.9, = 1.L and K = 1 and 7,(0) = 0.3, 
r9(0) = 1.5 and r3(0) = 0.5. Note that all densities are exponentially asymptotically 


stable. 
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Fig. 4. Trajectories and control function. 


In Figure. 5, %;,2) and zy the densities of tumor cells, predator hunting cells and 
resting predator cells respectively and equivalent control u,, in equilibrium point Fy, are 


134 M.R. Jafari, M.R. Zarrabi, S. Effati 


displayed against time for the system parameters and initial densities ag = 3, ag = 5, 
a3 = 4, a4 = 3, a5 = 2, ag = 5,01 = 0.9,a9 = 1.l and K = 0.15 and 271(0) = 0.3, 
r2(0) = 1.5 and r3(0) = 0.5. Note that all densities are exponentially asymptotically 
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6 Conclusion 


In this T the problem of instability of cancer self-remission and tumor system us- 

ing asin mode conto approach ave heen studied. The postive equim pon 

pap The stability and instability of the equilibrium points of 

are studied using the Lyapunov linearization approach. The eq srobont oobeaph fect for 
i of unstable equilibrium points ie tig and sliding surface are 

Analysis numerical examples for the controlled system were carried out for various pa- 

Tameters values and different initial densities. Matlab has been used for computations in 


this paper. 
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